Two loop expansion of the Schrodinger functional coupling a SF in SU (3) 
lattice gauge theory 

A. Bode a * 

a Institut fiir Physik, Humboldt-Universitat zu Berlin, 
Invalidenstr. 110, D-10115 Berlin, Germany 

The two loop coefficient of the expansion of the Schrodinger functional coupling in terms of the lattice coupling 
is calculated for the SU(3) Yang-Mills theory. This coefficient is required to relate lattice data to the MS-coupling. 
As a byproduct of the calculation, the Schrodinger functional is improved to two loop order and the three loop 
coefficient of the beta function in the SF-scheme is derived. 



1. Introduction 



2. Definition of a t 



In the framework of the Schrodinger functional 
the renormalized coupling a SF can be traced from 
low to high energy numerically on the lattice with 
finite size techniques. At high energy a conversion 
to etui's with perturbation theory is possible. The 
relation can be calculated by expanding both a-jjg- 
and a SF in go, the lattice bare coupling. 

This program is completed for SU(2) [||,f|,|| 
where the two loop relation between ojjs and a SF 
was required in order to avoid a significant source 
of error in the conversion to Ujjs- In the SU(3) 
case |(| the remaining step is to compute the two 
loop coefficient between a SF and go. 

This coefficient also enters into an extension of 
the program to include quenched quarks, which 
is discussed in Martin Liischers contribution to 
these proceedings. 

The calculation is analogous to the SU{2) case 
H). We have to calculate the perturbative coeffi- 
cients depending on the box size / = L/a, where 
a SF is defined. From this we are able to extract 
the continuum results and the O(a) improvement 
terms. Combining the numerical data and the 
two loop relations results in the expected error 
reduction comparable to the SU(2) case. 
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The coupling a SF {q) = g 2 SF {L) / {An), at q 
L^ 1 is defined via the effective action T: 



exp(-r) 



D[U] exp(-5[t/]) 



where the Wilson action 
9o „ 
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is modified by the weight w to achieve the O(a) 
improvement: w differs from one only for the pla- 
quettes attached to the boundary fields, where 
w = c t(<?o) is expanded in a series of g Q . The 
coupling 



(L) 



y'{l,t,) 



(3) 



r)=0 



is normalised via the classical action minimum 
Tq. The derivative is with respect to 77, which 
parametrises the boundary fields applied at t = 
0, L. The remaining three space dimensions have 
periodic boundary conditions. Details about the 
chosen diagonal and constant boundary fields can 
be found in 

3. Perturbative expansion and calculation 

The perturbative expansion has to be taken 
around the induced nonzero background field. 
Since the propagators are not known analytically, 
we calculate them numerically for each / and sum 



2 



them up within each of the 13 diagrams JT^]. The 
use of symmetries (for instance translation invari- 
ance and the cubic group in space) in this context 
reduces the numerical effort. The propagators are 
diagonal in the root basis of SU(3) and in spa- 
tial momentum space, and have been calculated 
there by a recurrence relation in time. Numerical 
efficiency and preserving maximal precision was 
taken into account in the program. In particular 
the 77 derivatives were taken analytically. 

The expected independence of the gauge pa- 
rameter Ao of T' 2 , the validity of the symmetries, 
a comparison of r' 2 with a numerical derivative of 
F 2 and the recalculation of the SU{2) results are 
applied as tests as well as an independent calcu- 
lation of Peter Weisz for small /. 

In 2 months CPU-time on a HP735 with 
128MB main memory we obtained m 2 (I) in the 
range 7 = 4... 32, where: 



mi (7) = ml(I) +c ( t 1) m b 1 (I) 
m 2 (I) = m a 2 {I) +c { pm b 2 {I) + 
+ c ( t 2) m 2 l (I)+m 1 (I) 



(4) 
(5) 

W) (6) 



From [|| ml is known and m\ , m 2 , m 2 are derived 
analytically. 

Symanzik's analysis suggests as asymptotic ex- 
pansion of m 2 



(/) = g rf +^ln(/)+^ln 2 (7) (?) 

n=0 



where we expect for instance t\ = from tree 
level improvement. 

Using the Liischer- Weisz || blocking method, 
which is based on 

R n [f](L) = f(L) + -^-f(L) n>l (8) 

With f(L) =a + b\n m (L)/L n we get: 

s f a + c/L n+1 m = 

R n [f](L) = \ ] m -i, n/Tn . n • (9) 
l^a + cln (L)/L n m > 

This enables us to cancel some powers in the 
residual terms. We used a symmetric lattice 



derivative and traced the numerical errors dur- 
ing the blocking. The error of the constant term 
was obtained by a fit to the residual terms. 

The universal /3-function coefficient b\ was sub- 
tracted after confirming it within 1.3%. At 
present all errors are very conservative. A more 
detailed error analysis will be presented in flCfi . 
The one loop improvement coefficient Cj was 
confirmed within errors. This should be consid- 
ered as an additional test of the calculation. 

Up to 0(I~ 2 ) times logarithms we get: 



mi(I) = 2&olnOO 
m 2 (J) = 2&i ln(J) 



.368282(11) 
.048085(63) 



mx{lf 



(10) 
(11) 



The continuum results are obtained by formally 
neglecting all negative powers of /, so eqn. 
(]l0|),(pd]) are the continuum results. 

4. Applying the expansion 

Together with the expansion of otjjs m 9o Q 
we get: 

ot-Ms{sq) = a SF +ci(s)a SF (q) 2 +c 2 (s)a SF (q) 3 (12) 
Cl (s) = -87t6 ln(s) + 1.25562(14) (13) 
c 2 (s) = ci(s) 2 - 32n 2 b 1 ln(s) + 1.197(10) (14) 

With bf^ @ we are able to quote b s 2 F = 
0.4827(88)/(4tt) 3 , where a fit to the data of 
numerical simulations gives b 2 s = 1.5(8)/(47r) 3 . 
There is no reason that the coefficient of the fit 
coincides with b 2 F , since a SF is traced nonpertur- 
bativly on the lattice. Nevertheless the order is 
the same and we expect as quoted in [|| the same 
error reduction as in the SU(2). 

The perturbative relation ( [l2|) between a SF 
and ajfs involves the scale factor s. Fixing them 
by demanding Ci(s) = or C2(s) to be minimal 
gives comparable results: 

a w (2.048g) = a SF (q) + 0.271(ll)a SF (g) 3 (15) 
a— (2.529g) = a SF (q) - 0.36895a SF (q) 2 (16) 
+ 0.135{ll)a SF (qf 

where the quantities quoted without error are in 
all digits significant. Note beside the scale factor 
the small coefficients. This perturbative relation 
between the two physical couplings should be safe 
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for small enough couplings. Using the smallest 
numerically determined coupling a SF we get for 
the SU(3) gauge theory: 

a— (14.5 GeV) = 0.1146(22)(2) (17) 
a— (78.2 GeV) = 0.08407(121)(5) (18) 

where the first error arises from the numerical un- 
certainties and the scale determination at low en- 
ergies. The second error results from the residual 
evolution and the conversion to the MS scheme. 
This should be compared with the results where 
C2 is unknown and &2 /f is used: 

a— (14.5 GeV) = 0.1145(23)(15) (19) 
ai7 s(78.2 GeV) = 0.08380(121) (59) (20) 

The error reduction of the conversion by a factor 
of 10 was the main goal of this calculation. 
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Evolution of the coupling a SF together with the 
perturbative evolution in 1 and 3 loop order. 



5. Relation to the bare coupling 

The relation of a SF to the bare coupling ao — 
<?o/(47r) involves a large scale factor and large co- 
efficients: 

a SF (14.06/a) = a + 4.178(ll)ag (21) 
a SF (17.36/a) = a -.36895c^ + 4.042(ll)ag (22) 



Using the tadpole improved coupling So = ao/P, 
we get a better behaved scale factor, but the co- 
efficients are still relative large. 

a SJ r(1.285/a) = a + 1.914(ll)ag (23) 
a SF (1.586/a) = a Q -.368955$ + 1.778(11)6$ (24) 

As an example we can insert &q = 0.12 (3 = 
6.5 = 6/s§,P = 0.6384) and get in eqn. (§J) a 
1-loop effect of 4.5% and a 2-loop effect of 2.6%. 
This casts doubt on the accuracy of tadpole im- 
provement in this region. 

6. Conclusion 

The present two loop lattice calculation com- 
pletes the connection between a SF and ctjjs in 
the SU(3) gauge theory. It has lead to a reduc- 
tion of systematic errors comparable to the case 
of £17(2). 
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